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We consider the Cauchy problem for a pseudo-differential operator which has a translation- 
invariant and analytic symbol. For a certain set of initial conditions, a formal solution is obtained 
by a perturbative expansion. The series so obtained can be re-expressed in terms of generalized 
Feynman graphs and Feynman rules. The logarithm of the solution can then be represented by 
a series containing only the connected Feynman graphs. Under some conditions, it is shown that 
the formal solution uniquely determines the real solution by means of Borel transforms. The 
formalism is then applied to probabilistic Levy distributions. Here, the Gaussian part of such 
a distribution is re-interpreted as a initial condition and a large diffusion expansion for Levy 
densities is obtained. It is outlined how this expansion can be used in statistical problems that 
involve Levy distributions. 
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1. Introduction and overview 

It is a well-known problem that there is no closed formula for the transition densities of 
general Levy processes. In fact, though there are many examples of Levy laws where the 
density is known (see, e.g., [1, 9, 19]), a general formula, even in the sense of series expan- 
sions, seems to be missing. Small-time expansions for pure jump Levy densities have been 
considered (e.g., in [2, 17]). But such expansions require the knowledge of convolution 
powers of the Levy measures which, in most cases, will be hard to calculate analytically. 
Expansions are available for stable laws (cf. [19, 21, 22] and references therein). 

It seems that there is no expansion formula for Levy densities in the case where the Levy 
process has moments of all orders that is simple enough to be applicable in statistics. 
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In this work, we provide such a formula in the sense of Borel summable asymptotic 
expansions in two cases. Either the Levy process has to be started in a initial distribution 
that is given by a density of special form or it can be started in a point but has to have 
a (large) diffusive part. The techniques employed stem from the renormalization group 
of statistical mechanics [18] and the recent proposal of a generalized Feynman graph 
calculus [5, 6]. Some of our ideas are similar to recent work of physicists motivated 
by Feynman path integration (cf. [9], Chapter 20.4.6, for a large diffusion expansion 
for one-dimensional and symmetric Levy densities). In [11], a different line of approach 
for obtaining asymptotics of densities given as solutions of pseudo-differential operators 
is presented. Furthermore, the expansion that we present is related, but not identical, 
to the well-known bootstrap expansion in statistics (cf., e.g., [3, 8] for an introduction 
to this nourishing field). However, it seems that non-symmetric densities on R d , Borel 
summability, the graphic representation and so-called linked cluster theorems have not 
been considered before. 

Most of our considerations are not necessarily restricted to Levy processes, but hold 
for fairly general convolution semigroups. We study the Cauchy problem of the partial 
(pscudo-) differential equation 

^ t {4>) = *(V)$tW, (t, <t>) e ]o, oo[ x R d , 

(i) 

$„(0 =f(<t>), 

for a given initial state / : R rf — > C and a (pseudo-) differential operator with constant 
coefficients ^(V). 

We consider the case where ^ is the symbol of a (pseudo-) differential operator, that 
is, has an expansion 

oo oo d 

*(£)=£< c(n) 'On = E E c&U.e* w 

n=0 n=0Xi,...,X„=l 

£= ,£d) € that converges on some neighborhood of and / is of the form 

/ = c- V , y(0)^E(A (p) ,O p = E E ^l.xjxi-ix, (3) 

P=0 p=0Xi,...,X p =l 

for 4> = ... , 4> d ) G K d , where p is even and (A® , 4> m ) P = V'\ v i ^xl-x^X, ■ ■ ■ 

4>x p > for 4> 7^ 0. We obtain a formal expansion of the solution of (1) in powers of 
A^ p ) that can be represented in terms of generalized Feynman graphs [5, 6]. Here, C^' 
iPxi—x ) ano - {*xl—x ) are symmetric tensors of pth degree over W 1 (their respec- 
tive components). (•, -) p is the canonical scalar product on (R d )® p . 
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For example, in the case where A*- p ) = unless p = 2 (this will be the most interesting 
case) and the symbol fulfills 1 = 0, the lowest orders in this expansion are 

i[00 +20+0 +2O 

(4) 

_l_ 4 (g,-^.^ _|_ 0_ —® 2 _|_ 4 ®-r<^> j + 3^ or der and higher in \( 2 > , 

where each graph stands for a numerical expression that can be calculated in terms of 
the basic parameters C^... x , A^... x , 4> an d t by an easy algorithm ("Feynman rule") 
(cf. the following section). 

This solution is the well-known graphical expansion for the renormalization group 
equation [18] for the case when ^(V) is a second-order partial differential operator (e.g., 
the Laplacian) and (1) is hence the heat equation. The recently discovered generalized 
Feynman graph calculus allows the generalization of this technique to a large class of 
symbols. This generalization is carried through in the present article, having in mind in 
particular Levy processes with generators given by symbols of the form 

*H0=i(a, £}-<£, + (e^-ljdrfc), (5) 

JR d \{0} 

which lead to equations (1) in the following jump-diffusion form (Levy-Ito decomposi- 
tion): 

f)& d 8 d B 2 

at d(f> x x^x 2 =i dcj)x 1 d(t>x 2 

+ z [ [* t (^ + V )-$ t (0)]dr(vj). (6) 

JK d \{0} 

Here, a € K d , D is a real and positive scmidefmite d x d matrix, z > and r is a probability 
measure on M. d \ {0} such that its Fourier transform is entire analytic. The well-known 
interpretation of a is as the drift vector. D determines the diffusion part, whereas z and 
r give the frequency and distribution of jumps, respectively. 

The restriction of the set of initial conditions in (3) is unwanted in many cases. It can 
be circumvented by splitting the symbol (5) into a pure jump and a pure diffusion part. 
After some time, the transition kernel of the heat semigroup generated by the diffusive 
part alone is of course a Gaussian density, which, up to normalization, is of the form given 
in (3). This simple observation makes it possible to develop a perturbativc calculus for 
densities of Levy processes starting in a point in the limit of large diffusion. Apart from 
this assumption, the perturbativc formalism developed in this article is fairly universal 

The assumption concerning is without loss of generality since a non-zero can be compen- 

sated by a shift of <f>. 
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and we therefore hope that it can contribute to a better handling of Levy distributions 
in statistics. 

The article is organized as follows. In Section 2, we develop a perturbative formalism 
for the solution of equation (1). In Section 2.1. we prove the differentiability of the 
solution $4 in the parameters that define the initial condition and conclude that 
the series so obtained is asymptotic. Section 2.2 gives an expansion of the perturbation 
series in the coefficients of VP which in the case of Levy processes are just the cumulants. 
The rearrangement in terms of generalized Feynman graphs is presented in 2.3. In the 
following Section 2.4, we obtain an expression for the logarithm of the solution <I>t in 
terms of connected Feynman graphs. This does not only reduce the number of Feynman 
graphs that one has to calculate, but is also of interest in connection with applications 
to the maximum likelihood method. A brief consideration of Borcl summability of the 
solution for the case that the logarithm of the initial condition, V , is quadratic concludes 
this section (Section 2.5). 

In Section 3, we apply the general results to the special case of Levy distributions. 
Section 3.1 presents the large diffusion expansion of Levy densities. In Section 3.2, we 
illustrate our method by a second-order calculation of the (non-normalized) Levy distri- 
butions using Pade resummation. The result is used for a calculation of densities and 
quantiles that we compare with the Monte Carlo simulation of a Levy process with 
compound or pure Poisson jump part. 

Some technical material that is presumably standard to some readers, but possibly not 
to others, is deferred to Appendices A.1-A.3. 

2. Generalized Feynman graphs and the convolution 
semigroup 

2.1. A solution in the sense of formal series 

After the previous rough description, let us now fix the mathematical details. Let the 
pseudo-differential operator '5 fulfill the following requirements: 

(PI) the associated symbol \&(— i£) is analytic for £ in a neighborhood of 0; 
(P2) 9?*(-i£) < clog(l+ || £ ||), where £ e R d and 3? denotes the real part; 
(P3) VaeNg, 3N a G N, b a > such that < b a (l + ||£|| ) N ° V£eR d . 

Here, D a h(<p) = g ^a d h((f>). We furthermore assume, without loss of generality, that 

\P(0) = 0. If this is not the case and ^(0) = b (i.e., if the associated Levy process undergoes 
killing or creation with rate 6), then the solution $ t of equation (1) can be obtained 
from the solution $ t of the equation (1) with "J replaced by W = \I r — ^(0) and setting 
$ t (^) = e t *(°) $*(</)). 

We take specific initial conditions that fulfill 

(11) m = *r v M- 

(12) V is a polynomial as given in equation (3). 
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Certainly, the restriction in the initial conditions by (II)— (12) is disturbing. Nevertheless, 
(II)— (12) are general enough to completely determine the convolution kernel v t , as one 
can construct approximating sequences of the Dirac delta distribution in that class (e.g., 
take V e (<p) = | 1| <p\\ 2 - f log(27te), e > 0). 

Let T : <S(R ) — > 5(R d ) be the Fourier transform on the space of Schwartz test functions 
over R d , that is, 

Hf)(0=/ m^dfp, feS(R d ). (7) 

JR d 

We denote the space of tempered distributions, that is, the topological dual of the 
Schwartz space S{R d ), by S'(R d ). By the duality given by the L 2 scalar product on 
R d with respect to the Lebesgue measure, S(R d ) is densely embedded in S'(R d ) and T 
extends uniquely to the latter space by continuity with respect to the weak topology [4] . 
We denote this extension by the same symbol T . 

By condition (P2), e***- -1 ^ is polynomially bounded for all t, that is, |e'*( -1 ^| < (1 + 
||£||) ct , hence e'*( -1 ^ is in 5'(M d ) for all i > and c > as in that condition. Furthermore, 
by conditions (II) and (12) - the positivity condition on A® in particular - fp = e~@ v 
is in the test function space S{R d ) for B> 0. Hence, for / = /i, the convolution $ t (</>) = 
v t * f(4>) = (vt, fc/>) is well defined, where v t G 6>'(R d ) is the inverse Fourier transform of 
e**'-*) a j S dU(<p) = f(<p-4>). 

Lemma 2.1. $t is the (unique) solution of the Cauchy problem (1). 

Proof. The proof is a standard argument based on the Fourier transform of tempered 
distributions. For the convenience of the reader, the details are given in Appendix A.l. □ 



Following Poincare, a formal power series X)m=o^ mam ^ s cane d an asymptotic series 
for the function h((3), j3 > 0, at /3 = if, for N e N, 



lim ft 

/3\o 



-N 



N 



(8) 



Here, we want to expand the solution (f> t in powers of V. In many cases, this expansion 
is not convergent. However, we show that such an expansion is asymptotic in the above 
sense. To this aim, let, for 3>0, <&f = vt * fp- Expansion in powers of V then means to 
expand in the auxiliary parameter f3 and then set 3 = 1. By definition, the series in V is 
asymptotic to <f>t if and only if the expansion exists and is asymptotic to $f . 
The moments of the distribution v t are defined as follows: 



(<PX! ■••<Px n )u t = 



d" 



neN, Xi, 



,x n e{i,...,d}. (9) 



For a polynomial P in tp 1 the expression (P((p)) v± is defined by linearity of the moments 
in the polynomial's coefficients. 
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Proposition 2.2. The V -expansion of <I> t is given by 



oo 



(-1) 



7 ) I 



rn ! 



V+(<p) = V(fp-it>). 



(10) 



m=0 



Equation (10) is to be understood in the sense that the right-hand side gives the asymptotic 
series for the left-hand side. 

Proof. The statement is evident in the sense of formal power series. That the right- 
hand side is asymptotic follows from the fact that $f (</>) is in C°°((0,oo)) with respect 
to (3 and for any m G N, $pz&t (0) can De continuously extended to j3 = such that 
lim^o 3pr$f (<P) = (-l) m (V^ m (<?))„ t (cf. Appendix A.l for the details). A posteriori, 

we conclude that the extension of $f is C°°([0, oo)) in j3 and apply Taylor's lemma to 
conclude the argument. □ 

In general, the asymptotic expansion is not convergent. This applies in particular to 
\&(— i£) given in (5), where ^> is the generator of a Levy process. We will return to this 
point in Section 2.5, where Borel summability is proved in an important case. By the 
fact that the expansion is asymptotic, it is already clear that the first TV terms in (10) 
give a good approximation to the actual solution for V small and </> not too large, an 
error estimate being provided by Taylor's formula. In order to extract reliable data for 
larger V, resummation techniques have to be applied. Some first steps in this direction 
are given in Section 3.2. 

Obviously, all partial differential operators of the form \&(V) = J2Z=o J2x ± x =i ^xl—x n x 
d vx 9 -d Vx of even order " with (-l) fi/2 (C (B) ,'0® s )n > for or of odd order n 
with C<") real and (-l)^-!)/ 2 ^" 1 ), V® fi_1 )n-i > for ^ fulfill the conditions 
(P1)-(P3). This, of course, includes the heat equation (n = 2). More interestingly, not 
only local generators can be treated, but also generators of Levy processes that have a 
symbol that is analytic at 0. 

Proposition 2.3. The generator (5) with r a probability measure with analytic Fourier 
transform at fulfills the conditions (P1)~(P3). 

Proof. (PI), analyticity at 0, is immediate from the conditions. The real part of (5) is 
obviously bounded from above by 2z, which establishes (P2). To see (P3), let us note that 
the partial derivatives of the third part in (5) arc of the form z J Rd ^ yip a e 1 ^' lp ' ) dr(tp). 
As r has analytic Fourier transform at zero, all moments of r exist. This implies differ- 
entiability everywhere and even that D a ^f(— i£) is uniformly bounded in £ £ M. d . □ 



For the sake of mathematical completeness, we also give a sufficient condition for the 
convergence of the expansion (10). Its practical value is limited, however, as it excludes 
the most interesting examples. 
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Proposition 2.4. If the symbol "J does not necessarily fulfill (P3) but fulfills (P3'): 
4 r (— i£) is entire analytic and 4 f (— i£) < K(l + ||S^|| max ), £ S C d , /or some A' > 0, then 
the expansion (10) is convergent. 

Proof. The condition implies that v t has an entire analytic Fourier transform |e**^| < 
e Kt Q tK HSVIImax^ Application of the Paley- Wiener-Schwartz theorem [14] implies that vt 
has support in a square of length 2tK centered at the origin and hence compact support. 
Let % be a test function with compact support and x|suppi/ t = 

1. Then $f (</>) = v t * 

fp((/>) = {vtiX ' e - ^ v *). As the series expansions in ft, /3 S C, of c~^ v * and all of its 
derivatives converge compactly on suppx, the mapping C 3 ft — > x • c _/3V * S 5(M d ) is 
entire analytic in ft. Since i/< S <S'(R d ), the same applies to $f (</>). □ 



2.2. Expansion in the operator's coefficients 

In this section, we define the truncated moments of the convolution kernel v t as partial 
derivatives of W - as the coefficients of "J" times t, in other words. A result known as 
the linked cluster theorem (sometimes also named the cumulant expansion) establishes a 
combinatorial connection between the moments {tpx t • • ■ L Px n )v t and the truncated mo- 
ments (tpxx ■ ■ 'i)3x„)I. The asymptotic expansion of the solution of the Cauchy problem 

(1) can thus be directly expressed in terms of the coefficients of the symbol C x \_, x , the 

time t and the coefficients A^.,, x of the logarithm of the initial condition. 

Definition 2.5. The truncated moments of the distribution v t are defined as 

(<p Xl ---vx n )Z t =t a " *(01e=o, neN, x 1 ,...,x n e{i,...,d}. (n) 

O 1 ?^! • • • ot,x n 

By (2), one obtains that the truncated moments are the coefficients of the (pseudo-) 
differential operator ^P, that is, 

(<p Xl ~-<Px n )Z=tC™.. Xn . (12) 

The following classical theorem gives the relation between the ordinary moments 
''-Px^vt and the truncated ones introduced above. 

Theorem 2.6 (Linked cluster). For n S N, we have 

ieP({i,.,n})!=i \jeik I v t 

/={/!,..., h} 



where, for a finite set A, V(A) stands for the set of all partitions I of A into non-empty 
disjoint subsets {I\, . . . , Ik}- 
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For the proof, we refer to the literature (e.g., [5, 16]). 

Proposition 2.7. The solution of the Cauchy problem 3>t(</)), <j) £ M. d , has the following 
asymptotic expansion: 

**(*) = E<=£- E E E 

™=o Pi,— ,Pm=o Kcn( Pl ,..., Pm ) xj,...,x* =i 

iev(n( Pl ,..., Pm )\K) 

e ^!...,x--a^;...,^ d4) 



x n (-^i)n^r 

(s,?)eJf (=i 

Here,ti(pi,...,p m )=\J™ =1 {(l,n),...,(p n ,ri)} andC h =cty b for I t = {(a u &i), ... , 

x 01 -x ag 

(a q ,b q )}. 

Proof. The proof is by a straightforward calculation using (11) and (12): 

' (-1)' 



**W = E H~ W»> 

z — ' m! r 



oo , , m p d d 

E l -1 i \ " \ ^ \ " Api) i(Pm) 

m=0 Pi,...,Pm=Ojfi,...,Xi 1 =l X™,...,X«^=1 

x ((fej - ^xj ) ' ' ' (Px^ - 4>x^ )■■■ (<Px™ - 4>x?) ■ ■ ■ (<£x™ m - 4>x™ m )) Vt 

E^- E E "• E <:!...,*. o») 

m=0 ' Pl,-,Pm=OXj,...,Xl 1 =l X™,...,Xm n =l 

>< e n (-<m( n ^r) 

XCn(pi,...,p m ) («,g)£JC \ (r,t)en(pi,...,Pm)\K 'ft 

00 / iim P d d 

E l-j) \" \ ^ \ " x (pi) iW 

m! 2^2^ 2^ A x 1 1 -xi l --- A xf-x Sri 

m=0 Pi,-,Pm=Ox i 1 ,...,Xi 1 =l x™,...x™ m =i 

>< e n maj) e n*^ 

xca(pi,...,p m ) {s,q)&K iev(Q(pi,...,p m )\K) i=i 

i={ii,...,i k } 
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□ 



Though the expansion is explicit, it is also rather messy. To simplify it, we introduce 
a graphical calculus in the next subsection. 

If d — 1, the expression (14) simplifies considerably; one can drop all subscripts X\, . . . 
and thereby obtain the asymptotic expansion 



A; 



oo , 1 s m p rn pH hPm / 

**<*) = E*=£- E n^ ;) E + 

m=0 pi,...,p m =0i=l k=0 

(16) 

piH hpm-fc <J 

E tq E Mii.-.yll 00,1 - 

9=1 h<l 2 <--<l q 8=1 
ilH hi<j=piH hp m — fc 



X 



Here, h q (li,... ,l q ) is a combinatorial factor given by the number of partitions of l\ + 
■ ■ ■ + l q objects in subsets with l±, I2, ■ . ., l q elements. It can be explicitly calculated as 
follows. Let {l s : s = 1, . . . , q} = {ii, . . . , i s } be the set of different values of the l\,...,l q 
and let q u = ${n = 1, . . . , q : l n = i u }, u = l,...,s. Then 

h q (h,...,l q )=( h + --- + l A^- n . (17) 

The expressions (14) and (16) simplify further if certain coefficients C^ n ' and vanish, 
for example, because of symmetries. As this is easier to understand in terms of graphs, 
we postpone the discussion of this point to the following subsection. 

In the case that / is a probability density and ^ is of the form (5), it is clear that ^t, as 
a convolution of a probability measure with a probability density, is a probability density 
for all t. Note that the first N terms of the expressions (14) or (16) give a polynomial in cf> 
and t. Obviously, this plays havoc with normalization conditions. We thus have to keep 
in mind that the approximation given by the first TV terms in the perturbation series 
is only good for t and <fi sufficiently small. In Sections 2.4, 3.1 and 3.2, we will discuss 
strategics to deal with this problem. 

On the other hand, the expansion in terms of polynomials in <p makes it simplest to 
integrate the probability density and to get a formula for the distribution function. Tak- 
ing, for example, the one-dimensional formula (16), we obtain, in the sense of asymptotic 
series, 

P(a <Z t <b) 

= £ IT" £ (" + "* + "-)Stt- 1181 

m=0 p!,...,p m =0/=l fc=o v 7 



Feynman graph representation of semigroups 



331 



PlH hp m -fc q 

9=1 h<l 2 <--<l q s=l 
JiH Kg=pH hPm— 

where Z 4 is a random variable with probability density function $ t , that is, a Levy process 
with generator \I/ and initial distribution /. Again, the first N terms of the right-hand 
side give a good approximation for a,b,t and \( p > sufficiently small. 



2.3. The Feynman graph representation 

In this subsection, we represent the solution of the Cauchy problem (14) as a sum over 
generalized Feynman graphs that are evaluated by a simple algorithm called "Feynman 
rules." In physics, the result is known as the Feynman graph representation for amputated 
Green's functions for the case when the generator is some kind of Laplacian, that is, 
^(^) = x 2 =i ^ x iAafalfe"' [18]- Combining this with the generalized Feynman 
graph calculus developed in [5, 6], we obtain the first main result of this article. 

A graph is a geometrical object which consists of vertices and edges, that is, lines 
that connect exactly two vertices. Here, we need graphs with non-directed edges that 
have vertices of different types. Vertices can be distinguishable or non-distinguishable, 
depending on their type, and they can have distinguishable or non-distinguishable legs. 
We use the term "leg" for the part of the edge meeting the vertex. All of these graph the- 
oretic notions are properly defined in Appendix A. 2. But for a first reading, an intuitive 
comprehension is sufficient. 

Our expansion is represented as a sum over generalized Feynman graphs. They have 
the following properties. 

Definition 2.8. A generalized amputated Feynman graph is a graph with three types of 
vertices, called inner full, 2 inner empty and outer empty vertices, respectively; see Table 
1. By definition full vertices are distinguishable and have distinguishable legs, whereas 
empty vertices are non- distinguishable and have non- distinguishable legs. 3 Outer empty 
are met by one edge only. Edges are non-directed and connect full and empty (inner and 
outer) vertices, but never connect two full or two empty vertices. 

Let m £ No . The set of generalized amputated Feynman graphs with m inner full ver- 
tices with pi,..., p m the number of legs of the inner full vertices such that pj < p and 
\(pj) ^ o, j = 1, . . . ,m, is denoted by F(m). 

The name generalized amputated Feynman graphs has been chosen to distinguish the 
graphs used here from classical Feynman graphs, but also from the generalized Feynman 
graphs with outer full vertices used in [5, 6] (this latter type gives the right combinatorics 
to calculate the moments of $ t ). To simplify the notions, we here use the term Feynman 

2 In physics this type is called an interaction vertex. 

3 Empty inner vertices can be distinguished from empty outer vertices, however. 
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graph when we mean a generalized amputated Feynman graph and there is no danger of 
confusion. 

Let m,pi, . . . ,p m G No- We construct a one to one correspondence between pairs (K; I), 
where K C f2(pi, . . . ,p m ) and I g V(£l(pi, . . . ,p m )\K), and the Feynman graphs F(m,pi, 
. . . ,p m ) with m inner full vertices with pi, . . . ,p m legs in the following way: 

• For j = 1, . . . , m, label the legs of the jth inner full vertex (1, j), . . . , (pj,j). In this 
way, we obviously obtain a one-to-one correspondence between the legs of the inner 
full vertices and the elements in f2(pi, . . . ,p m ). 

We first consider the case where the Feynman graph G is given and we construct (K, I) : 

• Define K C f2(pi, . . . ,p m ) as the set corresponding (under the above bijection) to 
those legs of inner full vertices in the Feynman graph G that run into an outer 
empty vertex. 

• Let k be the number of inner empty vertices in G. For q = 1, . . . , k, let I q be the 
subset in fi(pi, . . . ,p m ) corresponding to those legs of inner full vertices that run 
into the qth inner empty vertex. Finally, let I = {Ji, . . . ,Ik}- Then, evidently, I £ 

r(n( Pl ,..., Pm )\K). 

Conversely, let JfCSl(jj lr ..,j) m ) and I GjP(fi(pi, . . . ,p m ) \ K), I = {h, ■ ■ ■ ,h}, be 
given. We construct a Feynman graph G £ F(m,pi, . . . ,p m ): 

• Draw m inner full vertices with p\ , . . . , p m legs and k inner empty vertices with 
jj/i , . . . , jjifc legs and p = §K outer empty vertices. 

• Connect the legs of the qth inner empty vertex with the legs of the inner full vertices 
corresponding to the elements in I q . 

• Also connect the legs of outer empty vertices with the legs of the inner full vertices 
corresponding to the elements of K . The result is a Feynman graph G 6 F(m). 

For an illustration of the above construction, see Figure 1. We have established the 
following result. 

Lemma 2.9. Let m,pi, . . . ,p m £ No- There then exists a one-to-one correspondence be- 
tween {(K,I):K C fi(pi,...,p m ),J e V(£l(pi,...,p m ) \ K)} and F(m,pi, . . . ,p m ), the 
set of Feynman graphs with m inner full vertices with pi,.. . ,p m legs. 

The following algorithm in physics is named "Feynman rules" and gives a numerical 
value to each Feynman graph. 



Table 1. Different types of ver- 
tices 





Full 


Empty 


Inner 


• 


o 


Outer 


X 
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Definition 2.10. Let G £ F(m) be such that pi,.. . ,p m are the numbers of legs of the 
inner full vertices of G. Then, for (j> £ M. d and t £ (0, oo), V[G] = V[G](t,<j>) is a real 
number obtained by the following method 4 : 



1. 



to each leg of an inner full vertex, assign a value X a = Xj G {1, . . . , d}, with a the 
element in Q(pi, . . . ,p m ) corresponding to that leg; 

2. multiply each outer empty vertex with a factor (—4>x b ), where b € f2(pi, • ■ ■ ,Pm) is 
the element corresponding to a leg of the inner full vertex which is directly connected 
with that outer empty vertex; 

3. for each inner empty vertex with I legs, multiply with a factor tCj = tCx ai —x a > 
where J — {a\, . . . , a{\ C f2(pi, . . . ,p m ) is the subset corresponding to the I legs of 
inner full vertices that are directly connected to the inner empty vertex; 



(v ■ 

4. for the jth inner full vertex, j = 1, . . . , m, multiply with a factor Az.j 



X-. •••x rt 



ana sum 



up over J2 x i x 1 -r 




Figure 1. Construction of a generalized Feynman graph from the set K and the partition 

i={h,h,h}. 



4 We use the labeling of legs of inner full vertices by £7(pi, . . . ,p m ) introduced above, but any other 
labeling would do. 
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Combining Proposition 2.7, Lemma 2.9 and Definition 2.10, one thus gets the Feynman 
graph expansion of the convolution semigroup. 

Theorem 2.11. The solution of the Cauchy problem (1) has the following asymptotic 
expansion: 



Apart from the more efficient notation that comes with the use of Feynman graphs, 
Feynman graphs can be used to reduce the number of terms in the expansion. Note that 
the value V[G](t, 0) attributed to the graph G depends only on the topological Feynman 
graph. Hence, one can change the sum in (19) into a sum over topological Feynman graphs 
by multiplying V[G](t, <fi) with a weight factor or multiplicity, that is, with the number 
of Feynman graphs that give the same topological Feynman graph. Such multiplicities in 
low orders can be calculated by hand; see (4) for an example. 

Equation (19) is still fairly general. Often, a number of Feynman graphs do not give 
contributions. For example, in the case where = one has CW = for n odd 

and hence all Feynman graphs with at least one empty vertex with an odd number of 
legs can be omitted in the expansion. 

Another simplification occurs if the V((p) is quadratic in <p. Then, the only kind of 
inner full vertex that occurs is of the type -•- and can thus be considered as an edge 
of new type. Note that the legs of the inner full vertex are distinguishable, hence this 
new type of edge is directed, which changes multiplicities. Finally, if XjPx = ^XiX 2 is 
a multiple of the Kronecker symbol, the evaluation of a graph such as the one on the 
right-hand side of equation (4) can be simplified as follows. 

Corollary 2.12. IfV{4>) is quadratic, then the set of Feynman graphs F(m) of order 
m can be identified with the set of all graphs Q{m) with two types of indistinguishable 
vertices, called outer empty and inner empty, with to directed edges such that an outer 
empty vertex is hit by exactly one edge. 

If, in addition, A*- 2 - 1 is diagonal, 5 the value V[G](i, </>) attributed to such a graph G £ 
Q(m) can be calculated according to the following rules: 

1. to each of the m edges, assign a value X a £ {1, . . . , d}, a = 1, . . . , to; 

2. for each inner empty vertex with I legs, multiply with a factor tCx ai —x a > where 
a±,...,ai correspond to the numbers given to the edges that are connected to this 
vertex 6 ; 

3. for each edge, add a factor \, or, alternatively, add an overall factor X m ; 

4. sum up over all edges, that is, perform the sum 7 x =1 • 




(19) 



m=0 



G£F(m) 



5 If A( 2 ) is symmetric, this can always be achieved by choosing a suitable basis on R d . 

6 If an edge connects the vertex with itself, the corresponding index has to be repeated twice. 

7 Of course, for d = 1, the same simplifications as discussed at the end of Section 2.2 occur. 
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As a side remark, we note that the combinatorics of the Feynman graph formalism are 
dimension independent. Only Feynman rules depend on dimension. The number of steps 
needed in the evaluation of a given Feynman graph grow like ~ d m with the dimension, 
where to is the order of the graph. If one compares this with an evaluation of <3?t by 
a multidimensional Fourier transform, the steps needed for the evaluation grow rather 
like ~ m d , with to a parameter presenting, for example, the number of summands in an 
approximation of the one-dimensional Fourier integrals by a sum. For large dimensions, 
the Feynman graph formalism thus promises better results than numerical evaluation of 
Fourier transforms. Since the Feynman graph formalism has its roots in infinite dimen- 
sions, this should not be a surprise. 

2.4. The logarithm of the solution 

In this subsection, we obtain an asymptotic expansion of the logarithm of the solution 
in terms of connected Feynman graphs. In many points, the argument follows [5]. 

A graph G G F(m) is said to be connected if and only if each two vertices are connected 
by a walk passing through finitely many edges. We denote the collection of connected 
mth order Feynman graphs by F c (m). 

On the level of partitions, connectedness can be expressed as follows. 

Definition 2.13. Let m,pi, . . . ,p m G N be given and Jj = {(j, 1), . . . , (j,Pj)}, j = 
1, . . . , to. We say that the pair (K, I) with K C 0(pi, . . . ,p m ), I G V(fl(pi 1 . . . ,p m ) \ K), 
I = is connected if it fulfills the following condition: ^11 < i\, ...,i q < k, 

1 < q < k and 1 < j±, . . . , j s < m, 1 < s < m, such that 



Let V c (Q(pi, . . . ,p m ) \ K) be the collection of partitions / € V(fl(pi, . . . ,p m ) \ K) such 
that {K,I) is connected. 

Lemma 2.14. A graph G € F(m) is connected if and only if the corresponding pair 
(K, I) is connected. 

Proof. We only sketch the idea of the proof. In terms of the corresponding Feynman 
graph G, the condition in Definition 2.13 means that one cannot find a true subset Q\, Q2 
of empty and full vertices such that all edges from an empty vertex in Q\ end up at a 
full vertex in Qi and vice versa. This is evidently equivalent to the connectedness of G. 
For details, see [5]. □ 

Definition 2.15. Let Jj, j = 1, . . . , m, be as in Definition 2.13, (f> G M. d and X a = X\ G 
{1, . . . , d} for a = (j, I) G fi(pi, ■ ■ ■ ,p m ). We use the abbreviation Jj = Ylaejj (fx a — 4>x a ), 




(20) 



a=l 



a=l 
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j=l,...,m. The block truncated moment functions {Y]uLi jf)vt °f v t recursively 
(in m€N) defined as follows: 

k I \( T ) 

(4---Jtu= e n(n j ") • ( 2i ) 

IeV(l,...,m) 1=1 \ qeli I v t 

I={h,...,I k } 

The symbol (T) in Definition 2.15 means that each polynomial (in ip) j£ in the com- 
binatorics of truncation is treated as a single object. 

Proposition 2.16. Let J\, . . . , J m be as in Definition 2.13. We then obtain the following 
expansion of block truncated moments into ordinary truncated moments: 

{Ji---Ji)iV= E Ui-MllwZ- ( 22 ) 

ifcn(pi,..., Pm ) jei< i=i 

Proof. To verify equation (22), we have to insert the right-hand side of this equation 
into the defining equation for the left-hand side. We thus have to show that 



(Jt-j&u= e n 

IeV{l,...,m} 1=1 

I={h,...,I k } 



e iihmii^ 

Qi={Qi,u---,Qi,k,} 



(23) 



On the other hand, (jf ■ ■ ■ J^]v t can be directly expanded in terms of truncated mo- 
ments: 

k 

(4---Jtu= e nc-wnw- (24) 

Kcn( P1 ,..., Pm ) jeK i=i 
Rev(n( Pl .,..., Pm )\K) 

R={R u ...,R k } 

We have to prove the equality between the right-hand side of (23) and (24), that 
is, we have to construct a correspondence between the pairs (K, R) and the objects 
(J, (K±, Qi), . . . , (-K&, Qk)) indexing the sum in (23) such that their contributions to the 
sum are equal. 

Let (R,K) be given. We say that R connects q and j (in notation, q^R.j) if the full 
inner vertices corresponding to J q and Jj , respectively, are connected in the generalized 
Feynman graph corresponding to R. It is easy to see that R is an equivalence relation on 
{1, . . . , m}. Let / = {Ii, . . . ,Ik} be the equivalence class of ^_r. Then / e V({1, ■ ■ ■ , m}). 
We set Ki = (U, eJj J q f~\K), 1 = 1,..., k, and Qi = {S e R: S C U g€ii J g }- Using the 
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equivalence relation it is now easy to show that Qi <G 'P c ({J qeIi J q ); see [5] for the 
details. 

The converse construction is trivial: take K = ljf=i Kl an d R = Uf=i Qi- These two 
constructions give mappings between the index sets of the sums in (23) and (24) such 
that the summands arc equal. Furthermore, the mappings are inverses of each other, 
which establishes the one-to-one correspondence. □ 

We have now completed the preparations to prove the following theorem. 

Theorem 2.17. In the sense of asymptotic series, log$t(</>) has the following expansion 
in terms of connected Feynman graphs: 

log* t (0) = X;M- J2 V[G\(t,4>). (25) 
* — ' ml 

m=1 GSF c (m) 

Proof. Note that by the ordinary linked cluster theorem (Theorem 2.6), the following 
holds, in the sense of asymptotic series: 

°° ( — 1 ) m 

log* t (0) =log^,e- v *) = Yl [ - J r(VP% } - (26) 

ml v 

m—l 

The assertion thus follows by application of Proposition 2.16 and Lemma 2.14. □ 
The above result is useful in three respects. 

First, if is the symbol of a Levy process and the initial condition / = e~ v is a 
probability density, then the solution $t(4>) is a probability density for all t. We have 
seen in Proposition 2.7 and in equation (16) that the first N summands of the asymptotic 
expansion for $t(0) is a polynomial in <j) and thus non-normalizablc. However, we derived 
an asymptotic expansion for log$t(0) such that the first N terms are a polynomial 
Pn{4>) m 4> that goes to — oo for large 4> (this is often the case for N odd) and we obtain 
$ t (0) ^q^ p nW as the solution in iVth order perturbation theory. While the large <ft 
behavior of this approximate solution is obviously still problematic, 8 one has at least 
gained normalizability. 

Second, in statistical estimations by the maximum likelihood method, one needs an 
explicit formula for log$f (<^>), which is just the sum over connected Feynman graphs. 

Finally, the connectedness condition also reduces the number of terms in the expan- 
sion. For example, equation (4), only the 2nd, 3rd, 4th and 7th terms contribute to the 
logarithm of the solution. 

8 The probability that a random variable Zt distributed according <E>t takes large values <j> is suppressed 
much too strongly. To adjust the large <j> behavior, this first step has to be combined with rcsummation 
techniques (cf. Section 3.2). 
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2.5. Borel summability 

In certain cases, a function may be recovered from its asymptotic series expansion. This 
class of functions comprises the Borel summablc functions. A sufficient condition for a 
function to equal its Borel sum is the following version of Watson's theorem due to Sokal 
[20]. 

Theorem 2.18. Let h be analytic in the interior of the circle Cr := {£ £ C:3?C _1 > 
R^ 1 } (cf. Figure 2(a)) such that 



N-l 



h(0 - ]T a m C 



771=0 



<Ao"m\cr (27) 



uniformly in N and £ £ Cr, with constants A, a > 0. Under this assumption, one may 
show that B(t) := X)m=o a mT rn /m\ converges for |r| < 1/a and has an analytic contin- 
uation to the striplike region S a = {r £ C:dist(r,R + ) < 1/<t} (cf. Figure 2(b)) obeying 
the bound 

\B(r)\<KeM\r\/R) (28) 

uniformly in every S a i , with a' > a. Furthermore, h can be represented by the absolutely 
convergent integral 

1 f°° 

h(C) = 7 e~ T ^B(T)dr. (29) 
(, Jo 

Conversely, if B(t) is a function analytic in S a and there satisfying (28), then the func- 
tion defined by (29) is analytic in Cr, and satisfies (27) with a m = -^-^B(t) | t= o uniformly 
in every Cr/ with R' < R. 

If h is analytic on Coo = U_r>o Cr, = {( £ C : SRC > 0} and (27) holds on that domain, 
then of course the representation (29) also holds for £ £ C^ and for ( real and f3 = $t( £ 
M + in particular. The aim now is to apply this result to the series given in (19). It is 
easy to show analyticity, as seen in the following result. 

Lemma 2.19. The function 3 /3 — > $f (4>) ( cf. Section 2.1 ) has analytic continuation 
C 00 3C^$i(<t>), Vt£R + ,0£R d . 

Proof. Note that for C £ C^, f c = e^ v £ S(R d ) and hence $£(0) = v t * f c (<p) is well 
defined. To show complex differentiability in £, we now proceed as in the proof of 
Lemma A.l. □ 

In order to check (27), we need an estimate on the multiple ( derivatives of &i(4>) 
on Cr or even on C^. As &t(4>) is not known explicitly, such an estimate is difficult 
to obtain. But, in the case that 'J is the symbol of a Levy process, the probabilistic 
structure helps. 




Lemma 2.20. If ^ is the symbol of a Levy process and V(tp) > V<p £ K d ; then 

<{v^)„ t , CeCoo.teR+^eKf (30) 



jm 



dC 

Proof. We have \-^^{4>)\ = \(v t ,V™ -e _cv ^)| < (V™) Vt as v t is a probability measure 
and |e" cv *| <1. □ 

Lemma 2.20 allows us, in the given situation, to check (27) by Taylor's estimate for 
the remainder and the calculation of the .ZVth coefficient in the perturbation series. Being 
a sum over Feynman graphs, the coefficient can be dominated by the number of such 
Feynman graphs times the maximal value of V[G](i, cj>) over all these graphs. The latter is 
easy to control. We therefore start with an estimate on the number of Feynman graphs. 
Let us first assume that the shift <f> is zero. Wc furthermore assume that V is homogeneous 
of degree p. Then, only graphs without outer empty vertices and with inner full vertices 
with p legs contribute. As such graphs are in one-to-one correspondence with the number 
of partitions of pN objects, we have to control the number of such partitions. 

The number of partitions of m objects is known as mth Bell number and we denote it 
by b m . It has the asymptotic behavior (see, e.g., [12]) 

b m - m- 1 / 2 A(m) m+1 / 2 c A(m) -" i - 1 , (31) 

where A(m) is implicitly defined by A(m)lnA(m) = m. Using Stirling's formula AH ~ 
v / 27te' JV+1 / 2 - )lnJV-Ar , we get the following asymptotic formula for bpN /N\: 

— b- N 5_^-l/2 e (p-l)(iVlnAr-iV)-l/21niV-l e piV(lnp-lnlnA(piV)) e (l/2)lnA(piV)+A(pJV) ) ^ 
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from which we obtain that 

^b pN <A'K N N\ (33) 

is fulfilled for p < 2 and K sufficiently large. In the following, we restrict ourselves to 
p=2. 

Let us now consider a Feynman graph G £ F(N) without empty outer vertices. If the 
modulus of the coefficients Cj^ x is uniformly bounded by An n for A, k sufficiently 
large, wc obtain 

|V[G](i,0)| <max{A,l} N k N d N K 2N mnx{l,t} 2N VG £ F(N), (34) 

(2) 

with k — msocx 1 ,x 2 e{i....,d} I^XiX I* Thus, under the given conditions, the estimate (27) 
is fulfilled for = with h(() = $£(0), a m = { -^f- E Gef(m) V[G](t, 0), A — A' and a = 

dkn 2 K max{l, £} 2 max{A, 1}. It turns out that the case <fi ^ can be easily traced back 
to the case </> = 0. 

Theorem 2.21. Let ^ be a Levy symbol as in (5) such that supp(r) is compact and let 
f = c~ v withV(<p) = (\ ( - 2 \<p® 2 ) 2 . Then $<? (0), as a junction of C, fulfills the conditions 
of Theorem 2.18. In particular, the series 

B ( r >^)=Eti(2- E V[G](t,<t>) (35) 

™=0 ^ '' G<EF(m) 

converges on a ball of radius 1 jo centered at and has analytic continuation to the strip 
S a (for some a = a(t,4>) <oo). Denoting this continuation by the same symbol, we get 
the solution $ t of the Cauchy problem (1) via the formula 

poo 

®t{4>) = c- T B( T ,t,(p)&T. (36) 
Jo 

Proof. If supp(r) C {x £ M. d : \x\ < «}, then the uniform bound \C^... X \ < An n holds 

for A sufficiently large, since C Xi ... x f° r n > 2 up to a factor z is an nth moment of r. 
From the considerations above, the assertion now follows for <f> = 0. 

If <fi 7^ 0, wc modify for a fixed time t = to the Levy symbol by a deterministic 
term r^t^— i£) = i£) + i</> • £Ao generating the shift <j) at the time t . If we let 
©t,t o .0 be the solution of the Cauchy problem with ^ replaced by T to .^, then obviously 
^ > t ('/ > ) = ®t t <*(^) - Since the coefficients of \1> fulfill the above uniform bound, the same 
holds for the coefficients of Vt ,<j> (with A changed depending on to and </>). The assertion 
holds for 6j o tg ^(0), hence it also holds for $ to (0). 

In fact, the two expansions are equal. The only difference is that in the graphic ex- 
pansion of Ot ,t ,</)(0), there are no empty outer vertices. This is compensated for by 
evaluation of the empty inner vertices with one leg by (C x — <j>x)- Multiplying out the 
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parenthesis just generates the contributions of a non-modified inner empty vertex with 
one leg and the outer empty vertex. □ 

The restriction to p = 2 still covers a set of initial conditions / that determine the con- 
volution kernel vt completely (cf. Section 2.1). It is also sufficient for the most important 
application that we have in mind; see the following Section 3.1. 

The second technical restriction on the support of the jump distribution r is more se- 
vere. It can be overcome by considering generalized notions of Borcl summability (cf. [13]) 
that for the case of Levy processes can be adapted to a growth like (~ AH) 7 of the mo- 
ments of the jump distribution r. In applications, this matters if one wants to choose 
(36) or generalizations thereof as the point of departure of the numerical evaluation (cf. 
[10], Chapters 16.3-16.6) and one has to take care concerning the speed of convergence 
of the (generalized) Borcl transform. This program, as with all high precision numerics 
based on our proposal, is beyond the scope of the present article. 

It would be desirable to have a similar theorem for log $ t (</>) . The analyticity of 
log$j(</>) on some circle Cr is immediate. The bounds (27) for the perturbation sc- 
ries at zero is a corollary of the above considerations since the connectedness condition 
reduces the total number of Feynman graphs. In particular, the Borel sum B c (T,t,(j)) of 
log $((</)) has a non- vanishing radius of convergence. But, for the log-solution, the analog 
of Lemma 2.20 is missing and therefore the bounds on the series coefficients at ( = do 
not produce a uniform estimate (27) on some Cr, R > 0. So, the existence of an analytic 
continuation of B c (t, t, 4>) is not proved. Of course, one can repeatedly use the chain rule 
to obtain a bound on log$£ (0), but that generates an additional factor ~ on ■ So, more 
sophisticated estimates or more powerful versions of Theorem 2.18 are needed [13]. 

Here, we do not enter into this subject. Numerical procedures to evaluate the integral 
transform (29) are based on the convergence of the Borcl sum which permits resummation 
as, for example, the conformal mapping method [10], Chapter 16. Such methods, relying 
on the first N terms of the perturbation expansion only, automatically produce analytic 
continuations of an approximated B c (t, t, <f) on some striplike region that arc sufficiently 
bounded in order to render (29) absolutely convergent. 

3. Application to Levy distributions 
3.1. Large diffusion expansion of Levy densities 

It is more interesting to get explicit (approximate) formulae for the distribution 
of a Levy process starting at a fixed point (e.g., 0) than to determine the den- 
sity of the distribution for an initial condition / = c~ v with V as above. Though 
the densities of Levy distributions can be expressed in terms of more or less ele- 
mentary functions for many interesting examples (see, e.g., [1], Chapter 1), or can 
at least be expanded in a series in other cases [17], a generic formula is miss- 
ing. 



342 



H. Gottschalk, B. Smii and H. Thaler 



In this section, we present asymptotic expansions of (non-normalized) Levy densities 
and their logarithms, respectively, for the case where the generator of the process has 
a large diffusion constant fj, > 0. The formulae are generic in the sense that no detailed 
assumptions on the distribution of jumps, r, are required. That is, we consider the gen- 
erator (6) with DxiX 2 = ^x l x 2 i &x 1 ,x 2 being the Kronecker symbol, and we assume 
that /i is large and hence 1/fi is a small parameter in which we would like to expand. As 
the initial value problem for the density of the jump-diffusion process starting at zero, 
we get, for t,<j)& (0, oo) x R d , 

— *($ = /1 A* t ($- $>x-— / [* t fa + ¥O-*t(0]drfe>), 

ot £^ d<(>x Jr^\{o} 

(37) 

$ (<A) = <W), 



with So the Dirac measure with weight one at zero and A the Laplacian on M. d . The 
solution, <&j ', of the diffusive part alone is, of course, 

$f ff -( ( /,) = (47t/it)" rf / 2 e-^« 2 /( 4 ^= f^j ' e -^W\ (3=l/(4 fJ ,t),t>0. (38) 

Let v t = v 3 t ump be the inverse Fourier transform of e** J P ( _1 4), where nj/j um P is obtained 
from \& letting /.t = 0. Obviously, we obtain the solution of (37) by §t{4>) = ^t* $^ lff \ 
With V(<p) = ||0|| 2 , fp = e~ 0v one thus gets (§)- d / 2 $ t (0) = v t * f p {<j>). We note that the 
expansion of the right-hand side in powers of 1/ n is equivalent to the expansion in powers 
of (3 as the mth coefficient only differs by a factor l/(4i) m . Furthermore, from now on, 
we omit the multiplication with factors (3 in the evaluation rules for the Feynman graphs 
in order to obtain a more explicit (3 dependence. Combining this with Corollary 2.12 
and Theorem 2.17, we obtain the 1/ fi (large diffusion) expansion for the non-normalized 
(log-) density function. 

Theorem 3.1. In the sense of asymptotic series in (3, we obtain the following large 
diffusion expansion of the non-normalized solution (^) -d / 2 $ t , where 4>t solves (37): 



~ **(</>) =£^f- E V lG](t^). (39) 

™=0 GeQ(m) 



Furthermore, again in the sense of asymptotic series, we have 



log 



-d/2 



E^ E V[G]{tA), (40) 



m. 

m=1 GeQc{m) 



with Q c (m) the set of connected graphs in Q{m). 
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Combining of (39) with Theorem 2.21 yields the following convergent representation 
as a Borel transform. 

Corollary 3.2. Suppose that the support of the jump distribution r is compact. Let 
B{r,t,4>) be the analytic continuation of X)m=o fefp ^GeQ(m) MffK^ ^° some sfnp 
o/ the form S a . Then, 

poo 

$ 4 (</,)=7t- d / 2 /3 d / 2 ^ 1 / e- T l f3 B{ T) t^)&T. (41) 
Jo 

A comparison between (39)-(41) and other possible expansions known from statistical 
mechanics seems to be of interest. In particular, readers familiar with the theory of 
classical gases might find the convergent small z ("low activity") expansion [16] more 
suited than the small j3 ("high temperature") expansion. For a small z expansion, a 
graphical representation by Meyer graphs also exists. But there is a decisive disadvantage 
from the point of view of applications. The small [3 expansion produces expressions that 
directly depend on the quantities of statistical importance - z and the moments r„ 
of the jump distribution r. Hence, these quantities can be obtained from an empirical 
distribution by fits. This can be done to some extent nonparametrically, that is, without 
any a prior assumption on the jump distribution, as any sequence of moments of the 
jump measure depending on arbitrary parameter sets can be inserted into our formula. 
In the Meyer series, however, the jump distribution r enters via non-trivial r integrals 
over products of Meyer functions f(s, s', <j), <j)') = c~ ss ' — 1 which, except for a few 
cases with r particularly simple, can only be done numerically. 9 

Similar difficulties arise if one wants to solve the inverse Fourier integral of e**'- -1 ^. 
Note that this integral is of oscillatory type, so even its numerics is not really trivial. 

3.2. Second-order perturbation theory 

In this section, we give an illustration of how to actually calculate the 1/fx expansion of 
a Levy density and how to improve the result using resummation techniques. We take 
(40) as our starting point. Calculations are carried through up to second order in (3 only 
and we use the simplest possible resummation algorithm. Thus, the idea of this section 
is by no means to provide high-accuracy numerics, but rather to illustrate techniques 
up to second order of perturbation theory, that in higher orders will provide reasonable 
numerics. On the other hand, this does not mean that the second order calculation is 
not at all interesting from a numerical point of view. A more detailed discussion can 
be found below. Calculations to higher orders and more sophisticated resummation (see, 
e.g., [7, 9, 10]) are left for the future. 

In this subsection, we assume, without loss of generality, that ^j ump (— i^) has a van- 
ishing first coefficient If this is not the case, then this can be achieved by a shift 
<fi^ <fi — a (cf. the proof of Theorem 2.21). 

9 Note that the small-time expansion in [2, 17] is equivalent to the small z expansion. 
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Under these conditions, the second-order perturbation expansion in /3 for the left-hand 
side of (39) is given in (4). Extracting the connected graphs, we obtain 



log 



-d/2 



CO + 2< 



. 4 < g h *-^g> + 48h<>] + C(/3 3 ). 



(42) 



As in equation (4) , we omitted the evaluation symbol V and multiplied by the multiplic- 
ities of the topological graphs. 

Equation (42) should work for z and 1/fj, sufficiently small and 10 iwl, but the positive 
sign in front of the (3 2 term implies that the exponential of the second-order approxima- 
tion becomes non-normalizable (not LMntegrablc) if (3 (or z) passes a certain threshold. 
To extract reliable data for higher values of /?, resummation is required. 11 The idea 
of Pade approximants is to replace the Taylor polynomial by a rational function that 
has the same Taylor series as the polynomial. For definitions and details, we refer to 
Appendix A. 3. 

Following the standard calculation in A. 3, we obtain as the second-order Pade approx- 
imant 



log 



M<t>) 



1 + CO 



o 



+ (D([3 3 ), (43) 



<8 > « ® 



which, if exponentiated, gives a normalizablc approximation to $ t for t, z and (3 arbitrary. 
Evaluation for d=l with = zr n , r n = J R s n dr(s), yields an expression in terms of 
the time t, the activity (mean frequency of jumps in unit time) z and the 2nd, 3rd and 
4th moments of the jump distribution r: 



los 



71 



-d/2 



-(3(tzr 2 + (f> 2 ) 



-. f3 tzn+lt 2 z 2 r*+Atzr 2 q 
1 + 2 tzr 2 +4> 2 



■+Atzrz<j> 



(44) 



We have tested formula (44) by comparison with large samples of pseudo-random num- 
bers that simulate a compound Poisson process with additional diffusion, Z that is, 



Z = X(a, a 2 ) + siYi(zi) - s 2 Y 2 {z 2 ), 



(45) 



10 By a scaling of fi and z, one can always obtain t = 1. 

11 In fact, resummation of expansions in physics leads to results with impressive numerical precision 
(cf., e.g., [15], page 506). Often, validity of the resummed approximation is only undermined by phase 
transitions - which are, of course, absent on the discrete, finite space {1, . . . , d}, as considered here. 
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where X(a,a 2 ) has Gaussian distribution with mean a and standard deviation a = 
l/y/2j3. Yj(zj) is Poisson distributed with parameter Zj, j = 1,2. s\,s 2 give the length 
of positive/negative jumps, respectively. X, Yi,Y 2 are all independent. This choice is 
motivated by the fact that pseudo-random samples for Z can be obtained using stan- 
dard software routines. The parameters of the model are z = z± + z 2 and r = (z\S Sl + 
z 2 8- S2 ) / {z\ + z 2 ), and for the shift (/)—>■ (f> — a, one obtains a = z\S\ — z 2 s 2 . The nor- 
malization constant of the exponentiated right-hand side of (44) has been calculated 
numerically. The range of parameters has been 0.1 < /3 < 10 and 0.5 < z < 4, 1 < Sj < 10, 
J = 1,2. 

If vI>j um P(— i£) is symmetric under £ — > — £ (here, si = s 2 and z\ = Z2), the results of 
the second order calculation are very good, as long as jumps are not too large (< v^V)- 
However, the same is true for a best-fit Gaussian with mean a and variance zr 2 + 2/z. 

The situation changes if the distribution of Z becomes skew (e.g., by setting z 2 = s 2 = 
0). In this situation, more heavy tails are being formed. A typical result is displayed in 
Figure 3. 12 The overall precision of the second-order calculation is again comparable with 
a best-fit Gaussian model, but the second-order calculation typically predicts quantiles 
q a with a in the order of a few percent or above 95% with higher precision than the 
Gaussian model (cf. Figure 3(b)). As such quantiles are mostly used in statistics, the 
second order calculation already has some practical relevance. 

Appendix 

A.l. Proofs for Section 2.1 



Proof of Lemma 2.1. e"*^ 1 ^ -> 1 as t \ in the S'(R d ) topology, hence v t = T^ 1 x 
( e **(-i«)) jc-i(i) = S in S'(R d ) for t \ 0. Here, 8 is the Dirac measure with mass 
one at zero. Consequently, <£>((</») = Vt * /(</>) So* f{4>) = /(</>) as t \ 0. 

It remains to show that <& t is diffcrcntiable in t and fulfills equation (1). First, the 

differential quotient - — - — -jErp — - — — converges to \&(— i^e**^ -1 ^- 1 pointwise for all ^ g M. d . 
By condition (P2), we get 

< |*(-^)c** ( - i?) |(l + !|e!l) |t_ *' |c (46) 

and by (P2) and (P3), this implies that the differential quotient is uniformly polynomially 
bounded (for, say, \t — 1'\ < 1). It follows by dominated convergence that the differential 
quotient converges in S'(M. d ) as t' — ► t, which implies that T(yt) is differentiable in t in 

12 Actually, in this example, cr Q auss = \/2fl = 1.58 < 12 = 2 X jump length, hence the jump part for 
z\ = 2 is not small compared with the diffusive part, which provides a nice illustration of the power of 
resummation. 
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rel. frequency and densities 99 %-quantiles 




20 40 60 80 0.0 0.5 1.0 1.5 2.0 

<j> Activity z 

(a) (b) 

Figure 3. (a) Relative frequency of pseudo-random numbers (sample size n = 10 6 ), density 
of 2nd-order Pade approximant (44) (solid line) and a best Gaussian fit (dashed line) for 
z\ — 2, Z2 — 0, si = 6 and [3 — 0.2; (b) plot of 99%-sample quantiles of pseudo-random num- 
bers (triangles), predicted 99%-quantiles from (44) (circles) and a best Gaussian fit (+). Here, 
the data are as in (a), but z 1 varies from to 2. 



the S'(R d ) topology. By the continuity of T , the same applies to v t and hence $*(<?!>) is 
diffcrcntiable. 

It is immediate that F(-£Vt){0 = ^(-iOe**^^ , from which we get that J r (^$ t )(0 = 
, 3>(— i^)jF(<i>t)(^). From the fact that the time derivative can be taken out of the Fourier 
transform, uniqueness of the Fourier-transformed solution follows. The argument is com- 
pleted by application of the inverse Fourier transform. □ 
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Proof of Proposition 2.2. In the following, let ||/i|| a ,jv be a Schwartz norm where 
A G N, a G Nq, that is, 



\a,N = SUp (1 

4>GTS. d 



</>\\)"\D a h(<l>)\- 



(47) 



Lemma A.l. For [3 > 0, &t{4>) is infinitely differentiate in (3 and ^ Tr $f((/)) 
(-l) m vt*(V m f fi )(tt>). 



Proof. Since <&f ((f)) =Vt* fp(<f) an d v t € it sufhces to show that the mapping 



(0, oo) 3 0-/,, 



-pv 



€ 5(R d ) is to times diffcrentiable on (0, oo), to € N arbitrary, 



and g^r//3 = (— l)' m V m ff3. In fact, this statement is trivial for m = 0. Suppose that we 
have proven the statement up to to — 1 . We have to show that the differential quotient 
(-Ijm-iym-ij fe-fr ) converges to (-l)" 1 ^ as /3' -> /?. As the common prefactor 
i s a polynomial, it suffices to treat the case m = 1. For </> G R d , we have 



Hence, 



/3-/3' 



^(/g ~/r)d7 
/?-/?' 



(48) 



//3 - //?' 



/3-/3' 



,iV 



<A sup (1+11011)^- 



a'eNg:|a'|<|a 



<t>ei 



sup \D a ,f0((f>)-D a ,f T ((t>)\ 



(49) 



for A" and N a > sufficiently large. The derivative D a >fp((j)) is of the form P a i(4>, f3)fp((/)), 
where P a > (0, (3) is a polynomial in (3 and 0. Hence, for \f3 — (3'\ < 1, [3, (3' > e > 0, 

SUp |A*'//3(£)-A*'/tW| 
re(/3,/3') 



: SUP \P a ,(^f3)f^)-P a ,(^T)f T (4>)\ 



, (50) 



< sup \p a ,(^f3)-p a ,(^T)\f f} (4,)+ sup \P a i^T)UM- U4>)) \ 

< K a , |/3 - (3>\{1 + \\<t>\\f°> (fM) + \vme^ L e-^ P, ^ P H 

with L = mm.j.ezTxd V{cfj) - e{\(r\(j)®P) p , K a i , Q a i sufficiently large. Multiplying by (1 + 
11011)^=' and taking the supremum over <j> £ M. d implies that the right-hand side of (49) 
vanishes as (3' — > (3. This completes the proof. □ 



It remains to show that the derivatives of $f(</>) extend continuously to [3 — 0. We 
start with two technical lemmas. 



348 



H. Gottschalk, B. Smii and H. Thaler 



Let h r £ 5(M d ), r £ (0,oo), be a family of functions such that h r — ► h for r \ in 
5(R d ), where J Rd /id£ = 1. Further, let h e (£) = e^h^/e) and h$(£) be defined likewise. 
Let r(e) be a function such that r(e) — ► as e \ 0. 

Lemma A. 2. Let F :M. d — > C 6e a function that is polynomially bounded and has 
polynomially bounded partial derivatives. Let h e r and r(e) 6e defined as above. Then 
lim eN „/ Id FK)x^ E) K)d^F(0). 

Proof. We consider the following estimate: 



< 



F(eZ)K {£) (Odt;-F(0) 



(51) 



(F(e£)-F(0))h(O<lS 



The first term on the right-hand side can be estimated as follows: 



F(eZ)(h r(e) (Z)-h(Z))dZ 



< 



(i+iia) 



N 



d£\\K(e) ~ h\\ ,N- 



(52) 



Recall that F is polynomially bounded. Thus, for N sufficiently large, the integral on 
the right-hand side converges to |F(0)| J Rd (i + ||g||)iv d£ by dominated convergence. Hence, 
the right-hand side of (52) tends to zero as e \ 0. 

The second term on the right-hand side of (51), for < e < 1, can be estimated by 



(F(eO-F(0))h(Od£ 



JM d 



(VF(sO,OHO^ds 



(53) 



<eK (l + ||£||)XO|d£, 



for K and N sufficiently large, by the polynomial boundedness of the partial derivatives 
of F. Hence, the second term on the right-hand side of (51) also tends to zero as e \ 0. □ 

Lemma A.3. Let V e {ip) = (A® , ip® p ) p + YfpZl (Ae P) , ¥® P ) P such that A^ p) -»• as e \ 0, 

p = 0,...,p-l and f e = c- v - . Let f°(cp) = e^* '* 3 ®*^ . Then / £ /° m S(R d ) as 
e\0. 

Proof. We first assume a = 0, then |/° - f e \ < ^ \Qc~ sQ \dsf° < |Q|cl Q l/°, with Q a 
polynomial that has vanishing coefficients as e \ 0. Multiplying this by another polyno- 
mial and taking the supremum over M. d thus gives an expression that vanishes in that 
limit. Hence, ||/ £ - /°|| ,jv for e \ 0. 
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Next, let a £ Nq, a ^ 0. We note that each contribution to D a (fo — f e ) is of the form 
f e Q P or of the form S{f° - / e ). Q is of dc gree at most p—1 (potentially, Q — 0). P,Q 
and S are polynomials. The coefficients of P and Q are vanishing as e\0. Hence, we 
can conclude, as above, that ||/ e — f°\\ a ,N — > as e \ 0. □ 

We now conclude the proof of Proposition 2.2. By Planchcrcl's formula, (P2)-(P3) and 
Lemma A.l, we have 

jm r 

$f W = (_i ri , i *(^/ /3 ) W = (-i r / r^^Kjirfiv^di. (54) 

For a polynomial W%) = Ep= (^ (p) > V®*%> we set W/jfoO = £p=o P 1 - p/p (w(p),<p®r) p . 
Let = e -< A(ii) ^ S5i5 )p- an d let /0(<£) = e"^. By Lemma A.3, -> /° in S(R d ) as 

/3\0. Hence, ^(J?) -^^~Hf°) m 5 ( Rd ) and /r< ^ = /°(0) = 1. 

Furthermore, //^(p) = e -< 3V i>^) = e -{V^)^ /p <p) = f$(P 1/p (f). Inserting this into the 
right-hand side of (54), we obtain 

Jm r 

<&? (0) = (-i) m / r^^vfXf/^wev^Hfl de. (55) 



Note that by (P2)-(P3), ^(iV^e 4 **- -1 ^ as a function in £ is polynomially bounded with 
polynomially bounded derivatives. We can thus apply Lemma A. 2 and we get 

jm 

Hm— $? (0) = (-irV7(iV c )e**(-^| c=0 - (-l) m (^ m (^), t . (56) 
This completes the proof. □ 



A. 2. Notions from graph theory 

For the reader's convenience, we develop some graph-theoretic notions, some of them 
being non-standard, essentially following [6]. 

Given a set M , let M s 2 be the set of subsets of M containing two elements. A graph G 
over M is then a subset of M 2 . The points of M are called vertices of G and the elements 
(i.e., non-ordered pairs of M) {7711,7712} € G are called the edges (or bonds). Q{M) stands 
for the collection of all graphs over M. We could also replace M 2 with M 2 = M x M, 
the set of ordered pairs. We then obtain a graph with directed edges. 

Let N C M. The following construction explains the notion of a graph with distinguish- 
able legs (endpoints of edges) at the vertices in N. Let M be a set, 7T : M — > M a surjective 
map and tt 2 :Q(M) — > C/(M) the induced mapping on graph level defined by application 
of 7r on the endpoints of edges. Let G £ Q{M) be fixed and for m £ M, let p(m, G) be the 
number of edges connected with m. Let (Mcttg) be as above such that for m £ M \ N, 
tt7TQ 1 ({m}) = 1 and ^^({m}) =p(m,G) for m £ N . %A denotes the cardinality of the 
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set A. A leg m at a vertex m is a point in 7T G 1 ({m}). The class of graphs with distinguish- 
able legs at vertices in N associated with an ordinary graph G E Q(M) is then given by 
[G\N] = {G e (tt^-^G) : Vm G N, m G ^(m) ,${e G G : rh £ e} = 1}. The class of graphs 
with distinguishable legs at the vertices in N is then [£|iV](A/) = Ug65(m) 

Let L C M and s be in the permutation group of L, Perm(L). Obviously, s extends 
to a permutation on M by putting s(m) = m for m £ M \ L. Let s 2 :Q{M) — ► Q{M) be 
the induced mapping on graph level. We can then consider the equivalence classes under 
exchange of vertices in L, [Q:L](M) = Q(M) / (Peim(L)) 2 . Such an equivalence class is 
called a graph with indistinguishable vertices in L. It is easy to check that the above two 
constructions are compatible, that is, a vertex can belong to a set L of indistinguishable 
vertices with distinguishable or indistinguishable legs and a distinguishable vertex can 
have distinguishable or indistinguishable legs. 

A graph with several types of vertices has sets of vertices Mi, . . . , Mk and is a graph 
over M = {Ji=i M k- 

Finally, a topological graph is a graph with indistinguishable vertices in each class of 
vertices and with indistinguishable legs. The multiplicity of a topological graph, given 
some other class of graphs (generalized Feynman graphs in the present case), is the 
number of such graphs that coincide with the topological graph up to the labeling of 
indistinguishable vertices and legs. 

Let us also note that graph-theoretic notions like the connectedness of a graph with 
distinguishable legs G are meant in the sense that the associated "projected" graph G 
or even topological graph Gt op are connected. 13 



A. 3. Pade resummation 



The idea of Pade resummation is to replace the Taylor polynomial with a rational func- 
tion that has the given Taylor approximation at the origin but often has a better large 
(3 behavior. We essentially follow [10], Chapter 16.2. The [M/iV]-Pade approximant is 
defined as 

The M + N + 1 coefficients clq, . . . , am; b\ 1 . . . , 6/v arc determined by the equation 

M+N 

[M/N)(f3) = J2 h m r + 0{p N+M+1 ). (58) 

m=0 



G as a graph in Q{M) might be disconnected even if G is connected. 
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Multiplication with the denominator of [M/N] and comparison of coefficients leads to 
the system of equations 



bnhM-N+l + bN-lflM-N+2 + ■ ■ ■ + bofiM+i — 
oat^m + frw-i^M+i + • • • + boh]\i + N = 



&o = l 



and 



b h = g 
bih + b Q hi = a\ 



&o = l. 



b M ho + b M -ihi+ ■■■ +b h M = a M _ 
Solving for N — M — 1, as required in Section 3.2, one obtains 

hi 



bi = - 
a a = h 



hi' 



hnh 



h i 



a\ = h\ 

With ho,hi,Ii2 as in (42), this leads to 

lOO +2 



61 2 

a = 0, 

Gl = - [CZ^ + ®— ^Hg>] , 



(59) 



(60) 



(61) 



(62) 



from which (43) follows. 
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